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Abstract 

A new method has been developed to monitor the long-term 
calibration of the SBUV and TOMS instruments. It is based 
on the tact that the radiance in one channel can be expressed 
as a linear sum of the radiances in neighboring channels. 
Using simulated radiances for the SBUV and TOMS instru- 
ments, various scenarios of changes in instrument calibration 
are investigated. Results from sample processing of SBUV 
data are also presented. 

Introduction 

SBUV/TOMS instruments were launched aboard the 
Nimbus-7 spacecraft in October 1978. Since that time, the 
two instruments have been providing useful data for monitor- 
ing global columnar ozone. However, due to a hazardous 
space environment, one of the l§ey components of the instru- 
ment, a diffuser plate has been considerably degraded. This 
dill user plate allows a direct determination of the incident 
solar flux at the top of the atmosphere. Estimating the degra- 
dation is a major challenge in deriving the long term trend 
from the two instruments. If this degradation is not properly 
accounted for, the estimated ozone amount will be in error 
and a time dependent bias w'ould be apparent in the data. 
Four years ago, Bbartia and Taylor proposed a method to cor- 
rect the diffuser degradation using estimated ozone derived 
from different pairs of instrument wavelengths. This method 
requires that in the most favorable conditions (high sun and 
low- ozone in the equatorial region) the ozone derived from D 
pair (X3058A : X3125A) and A pair (A3 125 A : A3175A) be 
the same. Their method was further refined for the TOMS 
instrument by Herman et al. (1991) and was used for repro- 
cessing Version 6 of SBUV/TOMS data. Both methods 
assume constant values for diffuser error across the pairs. 

In this paper we describe another method to monitor the drift 
in the measured albedo due to the degradation of the diffuser 
plate. We define albedo as the ratio of the backscattered solar 
radiance in the direction 0 (polar angle) and <}> (azimuth 
angle) to the incident solar tlux F o in the direction 0 q at the 
top of the atmosphere. The method is based on the fact that 



Fig. 1 Normalized contribution functions for 8 of the 12 
wavelengths of the SBUV instrument. 
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the contribution functions for neighboring channels of the 
SBUV/TOMS instruments overlap one another This is 
shown graphically in Figure 1. The overlapping implies that 
the observed radiance in neighboring channels should be 
highly correlated. We have used this property to estimate the 
calibration drift in the SBUV/TOMS instruments. Unlike 
Herman et al.’s method which is applicable only for total 
ozone wavelengths, the method proposed here is very general 
and it is applicable to both the total ozone and the profile 
wavelengths. 

Method 

• The Equations 

If we assume that the contribution functions for channels j 
+1, j and j -1 overlap, then we can express the radiance (I) in 
channel j as a linear sum of radiances in channels j+1 and j-1 
as: 

- a o (0 g ,Rj>j) + 

where 0 o , R and Q are, respectively, the solar zenith angle, 
the surface reflectivity, and the columnar ozone amount in the 
atmosphere. Respectively, p and s, represent the surface pres- 
sure and the dependence on ozone vertical profile and a 0 and 
a are regression coefficients. 

Since the dynamic range of radiance over the SBUV wave- 
lengths spans over three orders of magnitude, the regression 
coefficients in Equation 1 would also show a large spread in 
their magnitudes. To reduce the large spread in the magni- 
tude, we used log of albedo instead of radiance (I) in the 
Equation 1. That is, 

AT/e.AQ'M - * Aj.i{6',Rj>j)* 

Nj.JB'.R.Qjjj) * *N n (B 0 ,R,Q,Pf) (2) 

where N is -lOOlogjQ (I/F) and A 0 , and A are new regression 
coefficients. (Simulation results showed that Equation 2 is an 
equally valid form of a linear relationship between the chan- 
nels). 

From Equation 2 it follows that, for a well calibrated instru- 
ment, the change in N value in channel j due to change in 
ozone amount is related to changes in N values in channels 
j+1 and j-1 by the following relation: 

AtyQ) - A J . i HN J .,(Q)+A J _ i *N H (Q) O) 

where for clarity we have omitted the other parameters from 
the arguments of N. 

• Estimation of Calibration Error 

To determine the instrument calibration error we assume that 
Nj t (the N value for channel j at time t can be written as: 

N P - AT r AA/«)-AAr/Q) ( 4 ) 


where ANj (€) and ANj (Q) respectively represent the change 
in N value due to changes in the instrument calibration and 
the ozone in the atmosphere; then from Equations (3), and (4) 
it follows that, 

* A n bN j'i ~ iN j* A j- 1 iN nV A j*i ( 5 ) 

where we have replaced ANj(c) by €j and 5Nj = (Nj - Nj t ) 
Equation (5) implies that if we know the error in channels j 
and j-I then we can determine the error in channel j+1. 
However, the accuracy of the estimate would depend on two 
things: a) how well the contribution functions overlap and b) 
the uncertainties in the estimate of errors in channels j and j- 1 

• Simulation of Radiances 

To test the concept described above, we simulated nadir radi- 
ance using Dave/Mateer’s radiative transler code (Dave 
1964). The code divides the model atmosphere into 101 lay- 
ers and solves the auxiliary equations of radiative transler by 
successive iteration. The code assumes a Lambertain surtace 
at the base of the atmosphere. We used mean SAGE pro! ties 
for equatorial region (10S-10N) from 1985 and computed 
nadir N values for all twelve channels of the SBUV instru- 
ment. The simulations were done for five values of solar 
zenith angle (18° (1) 22°), five values of reflectivity (0.22 
(0.01) 0.26), and six values "Of ozfine (250 (5) 275 DU). 
These values represent the nominal range of conditions over 
the equatorial region. 

We examined two scenarios of instrument calibration error. 
In the first scenario we modified the computed N values by 
subtracting a AN (e, X) value from each simulated N value. 
AN (€, X) was assumed to be a linear function of wavelength. 
The AN (€, X) values for longer wavelengths were very close 
to Herman et al. values for the SBUV instrument from 1987. 
In the second scenario we further subtracted a constant AN(c) 
(=0.5) from the N values of every other channel. We did this 
to test the stability and robustness of the algorithm lor an 
unusual pattern in the instrument behavior. We also exam- 
ined the impact of an unusual change in the profile shape on 
the estimate of the instrument error. For this purpose we first 
modified the SAGE profiles following the predictions of the 
photochemical models for mid-latitudes for a six year period 
from 1979 to 1985 (Watson et al., 1988) and then computed 
the nadir N values for the same range of solar zenith angle, 
reflectivity, and ozone amount. The computed N values were 
further modified by subtracting a AN (c, X) value from each 
simulated N value as in scenario one. The results ot these 
investigations and for the sample processing of SBUV 
Version 5 data are given below. 

Results 

Before we discuss the results for the different scenarios men- 
tioned above, it is instructive to examine the relationship 
between the N values for X3125A and X3175A (Figure 2). 
These wavelengths along with X3312A form the ‘A and B’ 
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Fig. 2 Relationship between ^3125 an ^ N 3175' The ozone 
amount in the model atmosphere was varied from 255 
DU to 275 DU. 


pair of wavelengths from which the total ozone is estimated 
Irom the SBUV/TOMS instruments. The N values were com- 
puled for solar zenith angle 0 Q = 22° and surface reflectivity 
R = 0.24 using the actual and the modified SAGE profiles 
trom 1985. Also shown are the N values with synthetic 
instrument calibration errors. The graphs in Figure 2 show 
that all N values for both the actual and the modified SAGE 
profiles (without instrument calibration error) fall on the 
same straight line implying that, for these wavelengths profile 
shape and ozone amount have no effect on the linear relation- 
ship between the two wavelengths. On the other hand, the N 


3 

1 

2 


o 

■a 

Jo 

3 


Wavelength (nra) 

Fig. 3(a) Estimate of calibration error for scenarios 1 and 2. 
The solid and dashed lines refer to the assumed input error. 
The filled circles and squares refer to the estimated error 
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values with synthetic instrument calibration error (Scenario 1) 
fall on a separate line suggesting that the instrument calibra- 
tion error can be estimated using the N value relationship 
between the channels. 

Figure 3(a) shows the retrieved error for scenario 1 (solid line 
with filled circles) and scenario 2 (dashed line with filled 
squares). The retrieved errors are practically identical to the 
assumed input error. The results for the third case ( i.e., esti- 
mating error in the presence of a significant profile shape) are 
shown in Figure 3(b). Here the initial instrument error is 
shown by a solid line and the final retrieved error by filled 
circles connected by the dashed line. In retrieving the instru- 


ment error, we have assumed that the errors for channels 1 1 
(X3312A) and 12 (^3398A) are known. For the 
SBUV/TOMS instrument they can be estimated by monitor- 
ing either the minimum reflectivity over the ocean or snow 
reflectivity over Greenland or Antarctica. Herman, et al. have 
shown that it can be estimated with an accuracy of 0.5 %. The 
results in Figure 3(b) suggest that we can predict instrument 
error reasonably well (with an error of less than 0.5%) up to 
X3125A and to about 1% lor A.3058 A. The large error for 
shorter wavelengths is due to sensitivity to profile shape. That 
is, the upper part of the modified profile is not representative 
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Fig. 3(b) Estimate of calibration error when ozone profile is 
significantly changed. The solid line refers to the true instru- 
ment error The filled circles ( squares ) represent the estimat- 
ed calibration error when the ozone profile shape is not (is) 
known. 


derived. The error is significantly reduced if we have a priori 
knowledge of the modified profile or the dependent data set 
includes all possible profiles. An example of a priori knowl- 
edge of the modified profile is also shown in Figure 3(b) 
(filled squares) which shows that the difference between the 
actual and the predicted error over all wavelengths is practi- 
cally zero. 


We have also applied this method to a sample of SBUV (V5) 
data. The sample was constructed by selecting one day of 
data from every month (around the 15th day) over the equator 
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Fig A Estimate of A pair calibration error as a function of 
time in Version 5 of SBUV archive data. 
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Conclusions 

In this paper we have described a simple method to estimate 
calibration drift in SBUV/TOMS instruments. Simulation 
results indicate that the method is robust and can accurately 
predict instrument error in TOMS instrument. Also with a 
priori knowledge of profile shape changes, the method can be 
used to monitor the drift in all channels of the SBUV instru- 
ment. An application of the method to the SBUV (Version 5) 
data yields slightly higher drift values for A pair ozone than 
Herman et al. However, we believe that the difference 
between the two values will decrease if we use a larger data 
base for computing the regression coefficients and correct the 
N values for small changes in reflectivity and solar zenith 
angle. Also, because of the simplicity and fewer assump- 
tions, the method can be easily used for determining calibra- 
tion drift in the SBUV/2 instruments on NOAA-9 and 
NOAA-1 1 and the TOMS instrument on the Meteor-3 space- 
craft. 
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